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(— ^ , Two important features of TeV-scale radiative seesaw models, in which tiny neutrino 

^\J ' masses are generated at the quantum level, are an extended scalar (Higgs) sector and 

the Majorana nature. We study phenomenological aspects of these models at the ILC. 
^ , It is found that the Majorana nature of the models can be tested directly via the 

electron-positron and electron-electron collision experiments. 

1 Introduction 

^ i' Current neutrino data give clear evidence for physics beyond the standard model (SM), 

JL , indicating that neutrinos have tiny masses of the 0.1 eV scale. Such tiny masses may be 

q_) ' generated from dimension-five effective operators ^v'^^j^v-'j^cfp (fp , where A is a mass scale, 

rl^ , Cij are dimensionless coefficients, and (ffi is the Higgs boson. The neutrino mass matrix 

can then be given by IvVj — Cij{(jP)'^ / K. As the vacuum expectation value (VEV) {(jP) 
of the Higgs boson is 0(100) GeV, the observed tiny neutrino masses are realized when 
^ ' (cij/h) ~ 0(10"^^) GeV~^. How can we naturally explain such a small number? In the 

>0 , tree-level seesaw scenario, the Majorana masses of right-handed neutrinos (= A) have to be 

[^ ' set much higher than the electroweak scale [1] . 

^— N , Quantum generation of neutrino masses can be an alternative to obtain (cy/A) ~ 

O(10~^'*) GeV~^. Thanks to the loop suppression factor, A in these models can be lower as 
compared to that in the tree-level seesaw models. Consequently, the tiny neutrino masses 

f^ ' would be explained in a natural way only by the TeV-scale dynamics. The original model 

along this line was proposed by Zee [2] , and some variations were considered O IH [Sj HI [7] , 
where neutrino masses are induced at the one-, two- or three- loop level. It must be charming 
that they are directly testable at the LHC and the ILC. 

A general feature in radiative seesaw models is the extended Higgs sector. Discovery 
of extra Higgs bosons and detailed measurements of their properties at future collider ex- 
periments can give partial evidence for the models. Previous works on phenomenology in 
radiative seesaw models O [9l [10] mainly discuss the constraint on the fiavor structure and 
coUider physics for the Higgs sector [TTl [T^J [131 [IH HI] • Next, the Majorana nature is another 
common feature. To induce tiny Majorana masses for neutrinos, we need lepton number 
violating interactions in the Higgs sector [2j [4] or right-handed neutrinos with TeV-scale 
Majorana masses [5l|6l[7]. When the future data would indicate an extended Higgs sector 
predicted by a specific radiative seesaw model, direct detection of the Majorana property at 
colliders must be a fatal probe to identify the model. 

In this talk, we discuss phenomenology of TeV-scale radiative seesaw models, in particular 
a possibility of detecting the Majorana nature at collider experiments [1^. We here discuss 
three typical radiative seesaw models as reference models; the model by Zee and Babu [4], 
that by Ma [6], and that in Ref. [7]. 
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Figure 1: Feynman diagrams for neutrino masses in the model by Zee-Babu [1] (left), that 
by Ma [5] (center) and that in Ref. [7] (right). 



2 Radiative seesaw models 

The Zee-Babu model 

In the model proposed in Rcf. 0j (we refer to as the Zee-Babu model), in addition 
to singly-charged singlet scalar bosons uj^, doubly-charged singlet scalar fields k^^ are 
introduced, both of which carry the lepton number of two unit. The neutrino mass matrix 
is generated at the two- loop level via the diagram in Fig. [1] (left). The universal scale 
of neutrino masses is determined by the two-loop suppression factor l/(167r^)^ and the 
lepton number violating parameter //. The charged lepton Yukawa coupling constants y^. 
(ye ^ JZ/j ^ J/t < 10^^) give an additional suppression factor. Thus, any of /y or gij in 
Fig. [1] (left) can be of 0(1) when niu: and nik are at the TeV scale. The flavor structure of 
the mass matrix is determined by the combination of the coupling constants fij and yigijyj- 
The flavor off-diagonal coupling constants fij and g.^ induce lepton flavor violation (LFV). 

In the scenario with hierarchical neutrino masses, fij satisfy /e^ — fer — fp.r/'^- The 
typical relative magnitudes among the coupling constants gtj can be ^^^ : 5^,- : grr — 1 '■ 
Tn^/mr : irn^lmT-Y' . For g^^ ~ 1, the neutrino data and the LFV data give the constraints 
such as TOfc > 770 GeV and m^j > 160 GeV [9]. On the other hand, the constraints on 
the couplings and masses are more stringent for the inverted neutrino mass hierarchy. The 
current data then gives m^ ~ 825 GeV for g^;^; ~ 1 [5] . One of the notable things in this 
case is the lower bound on sin^ 26*13, which is predicted as around 0.002 [8]. 

The Ma model 

The model in Rcf. [Bj, which we here refer to as the Ma model, is the simplest radiative 
seesaw model with right-handed neutrinos iV^, in which the discrete Z2 symmetry is intro- 
duced and its odd quantum number is assigned to N'^. The Higgs sector is composed of 
two Higgs doublet fields, one of which (S) is Z2 odd. As long as the Z2 symmetry is exact, 
the neutral components of S do not receive VEVs. We have one SM-like Higgs boson /i, 
and four physical Z2-odd scalar states; ^° (CP-even), ^,° (CP-odd) and ^^ as physical scalar 
states. This Z2 odd Higgs doublet is sometimes called as the inert Higgs doublet [TT] or the 
dark scalar doublet [12]. The LEP II limits are studied in this model in Ref. [TT] . 

The neutrino masses are generated at the one loop level via the diagram depicted in 
Fig. [1] (center), in which Z2 odd particles, ^^ and N'^, are in the loop. The universal scale 
for neutrino masses is determined by the one-loop suppression factor l/(167r^), the scalar 
coupling A5 and the mass Mjv° of the right-handed neutrinos. 
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In this model, there are two scenarios with respect to the DM candidate; i.e., the hghtest 
right-handed neutrino iV^ or the Hghtest Z2-odd neutral field (^° or ^^). For both cases, 
there are parameter regions where the neutrino data are adjustable without contradicting 
other phcnomcnological constraints [18]. In this talk, we consider the case where the dark 
doublet component ^° is the DM candidate. If the mass of the DM is around 50 GeV, the 
mass difference between ^° and £,^ is given by about 10 GeV when A5 ^ 10~^. In this case, 
hf ^ 10^^ (z = e,iJ.,T) are required to satisfy the neutrino data. The relic abundance of 
such DM is consistent with the WMAP data [13]. 

The AKS model 

In the model in Ref. [7], which we here refer to as the AKS model, it is intended that not 
only the tiny neutrino masses and DM but also baryon asymmetry of Universe are explained 
at the TeV scale. In addition to the TeV-scalc right-handed neutrinos N'^ {a = 1,2), the 
Higgs sector is composed of Z2-even two Higgs doublets $i (i = 1,2), where the physical 
component fields are H (CP-even), A (CP-odd), H^ and h (CP-even). and Z2-odd charged 
singlets S and a Z2-odd neutral real singlet r/^. The neutrino mass matrix is generated 
at the three-loop level via the diagram in Fig. [T] (right). The mass matrix has the three 
loop factor l/(167r^)'^ with additional suppression by the charged lepton Yukawa couplings yi. 
The electron associated couplings hl'^ and the scalar coupling k are of 0(1) for rrip^^ ^ 0(1) 
TeV. The Yukawa coupling constants hf are hierarchical as ft.e'^(— 0(1)) ^ /i^'^ ^ h].'^. 

The parameter sets which satisfy the current data from neutrino oscillation, LFV, relic 
abundances of DM and the condition for strongly first order electroweak phase transition 
are studied in Ref. [7]. To reproduce the neutrino data, the mass of H^ should be 100 - 200 
GeV. This is an important prediction of the model. In order to avoid the constraint from 
6 — ?► 57, the Yukawa interaction for the doublet fields takes the form of so-called Type-X [M] 
where only one of the doublets couples to leptons and the rest does to quarks [19]. The 
physics of the Type-X two Higgs doublet model (THDM) shows many distinctive features 
from the other type of extended Higgs sectors. For example, H and A decay mainly into 
T~^T~ when tan/3 > 3 and sin(/3 — a) ~ 1 [Ij. There are basically two DM candidates, 
77" and N^. We here consider the case where 77° is the DM [7]. Some characteristics of 77° 
are discussed in Refs. [20]. The coupling constant of S'^S~h is required to be of 0(1) for 
strongly first order phase transition, whose indirect effect appears in the quantum correction 
to the hhh coupling constant as a large deviation from the SM prediction [7]. As long as 
kinematically allowed, S^ decays via S^ -^ H^if by 100%. 

3 Phenomenology in radiative seesaw models at the LHC 

The existence of the extra Higgs bosons such as charged scalar bosons, which are a common 
feature of radiative seesaw models, can be tested at the LHC. Details of the properties of 
such extra Higgs bosons are strongly model dependent, so that we can distinguish models 
via detailed measurements of extra Higgs bosons. In addition, as the (SM-like) Higgs boson 
h is expected to be detected, its mass and decay properties are thoroughly measured |21| . 
The radiative seesaw models with a DM candidate can also be indirectly tested via the 
invisible decay of h as long as its branching ratio is more than about 25% for rrih = 120 
GeV with £ = 30 fb~^ [22]. The phcnomcnological analyses at the LHC in each model are 
in the htcrature [H dO] [El [H [15] . 

At the LHC, via the physics of extra scalar bosons such as (singly and/or doubly) charged 

LCWS/ILC2010 



Higgs bosons and CP-even Higgs bosons, the structure of the extended Higgs sector can be 
clarified to some extent. In addition, the invisible decay of the SM-like Higgs boson and the 
mass spectrum of the extra Higgs bosons would give important indication for a possibility 
to a radiative seesaw scenario. However, although they would be a strong indication of 
radiative seesaw models, one cannot conclude that such Higgs sector is of the radiative 
seesaw models. In order to further explore the possibility to such models, we have to explore 
the other common feature of radiative seesaw models, such as the Majorana nature. In the 
next section, we discuss a possibility of testing the Majorana nature at ILC experiments. 

4 Phenomenology in radiative seesaw models at the ILC 

At the ILC, properties of the Higgs sector can be measured with much better accuracy than 
at the LHC, so that we would be able to reconstruct the Higgs potential in any extended 
Higgs sector if kinematically accessible. Invisible decays of the Higgs boson can also be 
tested when the branching ratio B(h — > invisible) is larger than a few % |23j . Furthermore, 
the Majorana nature in radiative seesaw models; i.e., the existence of TeV scale right-handed 
Majorana neutrinos or that of LFV interaction, would also be tested. 

4.1 Electron-positron collisions 

In the pair production of charged scalar bosons at the e~^e~ collision, which appear in 
the radiative seesaw models (w+w~ in the Zee-Babu model, £,~^£,~ in the Ma model, and 
S'^S" (and H^H^) in the AKS model), there are diagrams of the t-channel exchange of 
left-handed neutrinos or right-handed neutrinos in addition to the usual Drell-Yan type s- 
channel diagrams. The contribution of these t-channcl diagrams is one of the discriminative 
features of radiative seesaw models, and no such contribution enters into the other extended 
Higgs models such as the THDM. These i-channel effects show specific dependences on the 
center-of-mass energy ^/s in proportion to log s in the production cross section, and enhances 
the production rates of the signal events for higher values of -y/s. The final states of produced 
charged scalar boson pairs are quite model dependent but with missing energy; 



e+e ^w+cj -^^le.^VLyL^ 


[The Zee-Babu model] 


e+e- ^ ec ^ w+'^*^w-'^*\%', ^ jjjjUjeL'^L)ee , 


[The Ma model] 


e+e- -^ S+S- -^ H+H-'n^rf -^ T+r-VLVLrfrf , 


[The AKS model] 



where underlined parts in the final states are observed as missing energy. 

The Ma model : The coupling constants h" (a = 1,2) are strongly constrained from neutrino 
data and LFV data. As a typical choice of parameters, we consider m^^ = 50 GeV, m^. = 60 
GeV, m^i - 100 GeV, m^yj, = "^iv| = 3 TeV, A5 = -1.8x 10"^ K,hf^, K ^ 10~^ in which 
the normal neutrino mass hierarchy is realized. Because ft," are very small for a TeV scale 
TTiATQ , the contribution of the t-channel diagrams to the signal e+e" — )■ £,+£,~ is much smaller 
than that from Drell-Yan type diagrams. For most of the possible values of /i" and m^a 
which satisfy the LFV and the neutrino data, the contribution of the t-channel diagrams is 
negligible. The production cross section of a charged Higgs pair S,+C~ is therefore similar 
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Figure 2: The jets invariant mass distributions of the production rates of the signal in 
the Ma model at ^/s = 500 GeV. left : The di-jct invariant mass M{jj) distribution of 



the signal e+e" — > S,^^^ — > jjfJ-'-'S.r^r fo'' "^5± = 100 GeV. right : M{jjjj) distribution of 
e+e~ -)> ^+^" -> W+W^^^S,^. -> jjjjC for m^± = 150 GeV. In addition to the rate from 
the signal process, those for main backgrounds are also shown. 



to that in the usual THDM: about 92 (10) fb for m^± = 100 (150) GcV at ^i == 500 GeV. 
The produced (^ decay into W^^*^S.ri- 

In Fig. 2 (left), we show the invariant mass distribution of the di-jet jj of the production 
cross section of the signal, e+e~ -^ ^+C" -^ W+*W-*^°^° -^ 3il^<%r for m^± = 100 GeV. 
The main backgrounds come from WW. The jj^i/i events from ZZ, 77, and Z7 can also 
be the backgrounds. A factor of 0.1 is multiplied to the rate of the jjii^i backgrounds for 
the miss-identification probability of a muon. The signal is significant around M{jj) ~30 
GeV. The invariant mass cut (such as 15 GeV< M{jj) < 40 GeV) is effective to reduce the 
backgrounds. For the numerical evaluation, we have used a package CalcHEP 2.5.4 |24j . 

For m^± > mw + w^^, on the other hand, the signal VF+W^~^J?^° can be measured by 
detecting the events of four jets with a missing energy. The main background comes from 
W^W~vv and ti. By the invariant mass cuts of two-jet pairs at the W boson mass, the 
biggest background from WW can be eliminated. In Fig. 2 {right), we show the invariant 
mass distribution of jjjj of the production cross sections of the signal and the backgrounds 
without any cut. A factor of 0.1 is multiplied to the rate of tt background, by which the 
probability of the lepton from a W that escapes from detection is approximately taken into 
account. The signal is already significant. The invariant mass cut {M{jjjj) < 300 GeV) 
gives an improvement for the signal/background ratio. 

The AKS model : We take a typical successful scenario for the neutrino data with the the 
normal mass hierarchy, the LFV data and the DM data as well as the condition for strongly 
first order phase transition [7]; m^ = 50 GeV, TO/f± =100 GeV, ms± = 400 GeV, to. 



-Nl 



TO 



m 



3 TeV, K - 0(1), hi 



2 » hi 



K » K 



hi, sin(/3 - a) = I, tan/3 = 10. 



Because h\'^ ^ C'(l), the contribution from the i-channel A''^ exchange diagrams to the 
production cross section of S^S^ dominate that from the Drell-Yan diagrams [7]. The cross 
section is about 87 fb for ?7i5± = 400 GeV at -y/i = 1 TeV. As the decay branching ratio of 



LCWS/ILC2010 



1000 



T T + missing energy] , — , 

e 




o 100 



600 800 

Vs [GeV] 



1200 




Figure 3: left : The cross sections of the signal, e+e^ -^ S^S^ -^ t^t^ (+ missing 
energy) , in the AKS model as a function of the collision energy ^/s. right : The differential 
cross section of the signal for ^/s = 1 TeV as a function of the angle of the direction of the 
outgoing T~ and the beam axis of incident electrons. In addition to the rate from the signal, 
those from backgrounds such as t^t^ , t^t^Vv and H^H^ are also shown. 

S^ -^ H^rf is 100% and that of H^ -)■ t^v is also almost 100% because of the Type-X 
THDM interaction for tan/3 = 10, the final state of the signal is T^T'Vvrfrp with almost 
the same rate as the parent S^S^ production. The main SM backgrounds are t'^t^ and 
T~^T~T'i'. The pair production of the doublet like charged Higgs boson H~^H~ can also 
be the background. As the signal rate dominantly comes from the f-channel diagram, it 
becomes larger for larger -y/s, while the main backgrounds except for TTVe^e arc smaller 
because they are dominantly s-channel processes (Fig. 3 (lejf)). At y/s = 1 TeV, the rate of 
the signal without cut is already large enough as compared to those of the backgrounds. It is 
expected that making appropriate kinematic cuts will improve the signal background ratio 
to a considerable extent. The y/s scan will help us to confirm that the signal rate comes 
from the t-channel diagrams. Fig. 3 (right) shows the differential cross section of the signal 
at y/s = 1 TeV as a function of cos 6^- , where 9^- is the angle between the direction of the 
outgoing T^ and the beam axis of incident electrons. The distribution of the background 
from TT is asymmetric, so that the angle cut for larger cos0^- reduces the backgrounds. 



4.2 Electron-electron collisions 

The ILC has a further advantage to test radiative seesaw models via the experiment at 
the e~e~ collision option, where dimension five operator of e~e~(/)"'"0"'", which is the sub- 
diagram of the loop diagrams for neutrino mass matrix. This direct test of the dimension 
five operator is essential to identify the radiative seesaw models. 

The Majorana nature in the Zee-Babu model is in the lepton number violating coupling 



constant /i of k^ 



' , which generates the dimension five operator of e 



,+,.,+ 



tree level via the s-channel k exchange. The cross section of e e 



a{e e — > w w 
where F;, is the total width of fc^^ 
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is given by 
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s (s — m^)2 + m^F^. ' 
In the Ma model and the AKS model the operator 
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comes from the i-channel N"^ exchange diagram. The cross section is evaluated as 



(T(e e 
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■m^o 



t - m2, 



"N^ 
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1000 
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where n is the number of generation of right-handed neutrinos, (j)~ represents the Z2-odd 
charged scalar boson ^~ in the Ma model and 5*^ in the AKS model. The constants c" 
represent /i" or /i" in the Ma model or the AKS model, respectively. Due to the Majorana 
nature of the t-channel diagram, we obtain much larger cross section in the e~ e~ collision 
than in the e^e~ collision in each model. 

The mass matrix of left-handed neutrinos 
is generated at the one, two and three loop 
levels in the Ma model, the Zee-Babu model 
and the AKS model, respectively. Therefore, 
the coupling constants can be basically hierar- 
chical among the models, so are the cross sec- 
tions. For the typical scenarios in these mod- 
els, the cross sections are shown in Fig. 4. The 
rate of lo^uj^ in the Zee-Babu model can be 
larger than several times 100 fb for 800 GeV 
< yi < 1.5 TeV, where m^ = 1200 GeV, 
jji = 800 GeV and gee = 0.17 are asuumed, 
and Ffe is computed as about 168 GeV in our 
scenario (see [l^). It becomes maximal (sev- 
eral times pb) at -^/s ^ m/j, and above that 
asymptotically reduces by 1/s. The maximal 
value of the cross section is sensitive to the 
value of Qee and /i. The signal should be like- 
sign dilepton with a missing energy. On the 
other hand, in the Ma model, production cross sections of e^e^ — > S,^ ^'^ are smaller than 
10^'^ fb. Hence, most of the successful scenarios in the Ma model the process e^e^ -^ £,£,^ 
is difficult to be seen. In the AKS model, the cross section of e~e~ — > S~S~ is large, and its 
value amounts to about 15 pb at ^/s = 1 TeV in the present scenario. Above the threshold, 
the magnitude of the cross sections are not sensitive to ^/s, so that even if mg± would be 
at the TeV scale, we might be able to test it at future multi-TeV linear colliders, such as 
the Compact Linear Collider [25]. Because B{S^ -^ rfH^) ~ B{H^ -^ t^v) ~ 100 %, 
the signal should be t^t^vvi^i^ with almost the same rate as long as "t.o± < mN°- The 
background mainly comes from W~W~VeVe^ and the cross section is about 2.3 fb (22 fb) 
for Y^ = 500 GeV (1 TeV). The branching ratio for the leptonic decay of W bosons is 30%, 
so that the rate of the final state U'vvvv is at most 2 fb. Therefore, the signal in the AKS 
model and in the Zee-Babu model can be seen. 

There are several other models with lepton number violating interactions or right-handed 
Majorana neutrinos. Atwood et al. have discussed the signature of heavy Majorana neu- 
trinos in the model without Zi symmetry via charged Higgs pair production at e"'"e~ and 
e~ e~ collisions |26| . In supersymmetric models, Majorana particles also appear, and their 
effects also give similar i-channel contributions in the slepton pair production e^ er — > e^e~, 



Figure 4: Cross sections of like-sign charged 
Higgs pair productions in the Zee-Babu 
model and in the AKS model. 



whose cross section is at most 0(100) fb. The final state would be e 



e XX 



"^ for example. 



LCWS/ILC2010 



5 Conclusion 

We have discussed general features of TcV-scale radiative seesaw models. They are charac- 
terized by an extended scalar (Higgs) sector and the Majorana nature. Various phenomeno- 
logical aspects especially in experiments at the ILC. have been discussed in the Zee-Babu 
model (2-loop), the Ma model (1-loop), and the 3-loop model in Ref. [7]. While the ex- 
tended Higgs sector can be explored at the LHC, the Majorana nature of the models can 
directly be tested at the ILC via the pair production of the charged scalar bosons at the 
electron-positron and electron-electron collision experiments. 
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